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It was found recently that tunneling probabilities over a barrier is roughly 
twice as large as that given by standard WKB formula. Here we explained 
how this come from and showed that WKB method does give a good approx- 
imation over almost entire energy range provided that we use appropriate 
connection relations. 



PACS numbers: 02.60+y, 03.65Sq, 02.70+d. 
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I. Introduction 

WKB method was first invented by Jeffreys 1 and was applied to solve Schrodinger 
equation by Wentzel, Kramers and Brillouin 2 . It is a powerful tool and has 
many applications, e. g. , waves in a inhomogeneous plasma 3 . Consider a 
second order ordinary differential equation: 



where prime denotes differentiation with respect to argument. The function 
k 2 may be negative for some x. In this case, we define k 2 (x) = —k 2 (x). WKB 
method approximates solutions of (1) by combinations of 



is satisfied. If (3) is not satisfied over a range (or at a point), then we need 
to connect WKB solutions of both sides appropriately by some connection 
relations. One common situation where (3) breaks down is that k 2 near 
a turning point. 

In a recent paper 4 , a quantum mechanical tunneling example is consid- 
ered. The Schrodinger equation in this one dimensional problem can be 



i)"{x) + k 2 {x)i){x) = 0, 



(1) 





provided that 



k'/k 2 < I 



(3) 
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represented by (1) with 



k 2 ( , [E-a* forO<*<6 

y ' \ E otherwise, 1 ' 



where b = ^ E^/a and E, E\> and a are constants. E may be assigned 
the physical meaning of energy. It was shown that, for E < E b , the true 
(numerical) tunneling probabilities T n (see Appendix for numerical method) 
are roughly twice as large as the well known WKB tunneling factor: 

T w = e - 2 £ = ^{^hV^I-zV^^ (5) 

where a = \J E/a and j3 = \J ' E h / ' E. Values of T w and T n for different E are 
shown in Table I and Table II. 

This result was unexpected to us at first. After a little survey over some 
standard textbooks, it was found that two books did warn about such an 
error 5 ' 6 , although without elaborations. Another book did an example with 
correct connection relation 7 , but did not give general discussions. Two other 
books used the standard tunneling factor without warning 8,9 , although one 
of them indicated that it is just a rough approximation 9 . In this case, we 
have reasons to assume that the restrictions in applying the standard WKB 
formula are not very well known to students. Therefore, we would like to 
show how to apply WKB method to solve a tunneling problem more carefully. 

In section II, we briefly review the derivation of the standard WKB tun- 
neling factor using standard connection relations, which are also useful for 



our discussions later, over a typical turning point. In section III, the tun- 
neling factor is corrected when one turning point is a sharp edge. In section 
IV, further correction is made when the tunneling probability is not small. 
Section V considers the case that E — > 0. Section VI deals with the case 
that E E b . Finally section VII treats the case that E > E b . It will be 
shown that WKB method does give good approximation to the tunneling 
probability if we use different connection relations for different energy range 
appropriately. 

II. Standard WKB tunneling formula 

Consider (1) with x near a turning point x — 0. Assume that k 2 (x) can be 
approximated by: 

(6) 



k 2 (x) w ^- 
dx 



x=0 

Let (dk 2 /dx) Q > for the time being. Solutions of (1) with k 2 given by (6) 
can be written as combinations of Airy's functions 10 Ai(— Xx) and Bi(— Xx) 
with A = \(dk 2 /dx)^ 1 ^ 3 . Then the asymptotic behaviors of Ai and Bi give 
the standard WKB connection relations 10 : 

— OO i X > OO 

^e-/>* <- yfAi(-Ax) j-sin(JSkdy+z) (?) 

^J> dv JlBi(-Xx) — ). ^cos(/**di/+f). 
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For (dk 2 /dx) < 0, we only need to change A to —A in (7), interchange 
the limits of integrations and interchange left and right hand sides. These 
standard connection relations are valid only if the turning point is smooth 
enough, i.e., we require (6) to be valid up to some x where (3) is also satisfied. 
This is not always possible 3 . 

Consider a barrier with k 2 > for x < a and x > b, k 2 < for a < x < b 
and that WKB condition (3) is valid over the whole range of x except near the 
two smooth turning points a and b. Then for a incident wave from x = — oo 
direction, we have only outgoing wave for x > b: 



^ = ±=jS: kd y. (s) 



By (7), ip connects to: 



JL e f b x <^y + * e -f b x 



e~ in/A (9) 



^£ = e _ />*- il/4 ) (10) 



for a < x < b. Then for x close to a, 

J 

where I = k dy and is assumed large. By (7) again, ip connects to 

ip = 2 ^_ sin (J kdy + Ti/Aj (11) 

for x < a. Since the sine term is just a combination of incoming and reflected 
wave with equal intensity, we have reflection coefficient R w ~ 1 and tunneling 
probability T w = e~ 21 . This is how (5) come from. 



From the above derivation, we see that in order to apply the standard 
WKB tunneling factor (5), three conditions must be satisfied : 

(i) there are two and only two turning points; 

(ii) (3) is valid except near the two turning points where k 2 can be approx- 
imated by (6); and 

(iii) T w < 1 . 

We immediately see that the potential given by (4) does not satisfy (ii) since 
k 2 is discontinuous at x — b. This is the main reason for the factor of two 
error. Also, for larger E, T w is not very small, (iii) is also violated. 

III. Sharp turning point correction 

Let us now consider k 2 with a discontinuity at x — b, so that k 2 = —k 2 < 
at b — + and k 2 = k 2 > at b + + . We still assume that the WKB condition 
(3) remains valid except at b and near the smooth turning point x — a. For 
x > b, there is only outgoing wave (8). It connects to 




for a < x < b. We can check whether (12) is correct by substituting x = b 
into (8) and (12). We should see that ip and ip' are indeed continuous. Similar 
to steps (9) to (11) and still assuming I to be large, the corrected tunneling 
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probability can be found as: 



T i - K 2 I p T «- ( 13 ) 



For k 2 given by (4), 



AJ E b /E - 1 

T > = J 4ib- T -- (14 > 



It can be easily shown by (13) that T x < 2T W . The equal sign holds for 
kb = Kb. This explains the factor of two found by ref. 4, Some values of 
T w and T\ are shown in Table I to compare with true (numerical) tunneling 
probabilities T n . Parameters a = 0.040965, E b = 1.2776 were chosen so that 
we may compare with Table I of ref. 4. We see that T x is much closer to T n 
from E ps 0.2 to E ps 0.6 while T w is nearly a factor of two smaller. However, 
for large E when T w lager than 0.2, 7\ also fails since we assumed T w to be 
small in the above derivation. This error will be corrected in next section. 

IV. Finite T w correction 

First, let us consider the correction of T w itself. If we keep both terms in (9), 
then ip will connects to 

1 



2 sin(0) + l -T w cos(#) 



e~ in/4 (15) 



for x < a where 9 = /" kdy+n/A. Then by finding the coefficient of incoming 
term, the corrected tunneling probability can be found as: 

T w2 = T w /(l + T w /4) 2 . (16) 



Similarly, keeping both terms in (12) gives correction of T\ : 

T 2 = 71/ (1 + ^/2 + 7^/4). (17) 

Some values of T 2 are also shown in Table I and Table II. We see that T 2 — > T\ 
as T w — > 0. But for large T W) T 2 differs from T\ quite a lot and gives better 
approximation to T n up to E 1, or T w 0.65, which is quite large. 
However, for 75 — >■ T 2 also fails. For example, for E = E b , both T\ and 
T 2 equal to zero while T n is actually quite large. The reason for this error is 
that we assumed that WKB condition (3) is satisfied up to x = b. This is 
not true for E — > E^, since k is small, even for x = b, while k' is not small. 
So we cannot use (13) or (17) for E close to Eb. Another way to apply WKB 
method for this case is discussed in section VI. 

In the other extreme, the fact that for E — > 0, T w — ^-finite value while 
T n y 0, is another limitation of the standard WKB tunneling factor 4 . Al- 
though Ti and T 2 give the right value, i.e. zero, at E = as shown in Table 
I, the dependence of E is wrong as we can see from Table II. We see that 
T n oc E as E — > while T 2 m Ti oc \/E. So, there are extremely large error 
between T n with Ti or T 2 . We will discuss this case in the next section. 
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V. For E -> 



Consider the general case first. Let k 2 = E — > for x < a. By (12) and 
connection formula (7), we know 

1 



h V K b 



Ai(Ax) 



(18) 



for x — > a + , where A 3 = — (dk 2 / dx) x=a . Again, T w small was assumed, 
although this assumption can be removed if higher accuracy is desired. Now 
assume 

rj} = Ce iV ^ x + De~ iV ^ x (19) 

for x < a — > 0, where C and D are constants. Matching (18) and (19) at 
x = and using the fact that E — > 0, we got the tunneling probability: 

i Ve 



T 



>7i. 



|C|V£ 7rA [Ai'(0)] z 
Note that power series expansions for Airy's functions are 10 : 
Ai(z) = Cif(z) - c 2 g(z), 
Bi(z) = V3[ Cl f(z) + c 2 g(z)}, 

where 



(20) 



(21) 



/(*) 




z 3 

5-6 


1 + 


z 3 

8-9 


M 




\ z " 


1 + 


Cl 


= 0.355028, 








c 2 


= 0.258819. 









9-10 



So, Ai'(O) = — C2. In the derivation of (20), we assumed A finite and k 2 
smooth enough near x = a. However, for k given by (4), A = (AaE) 1 / 6 ->■ 0. 
This means that T oc E 5 ^ 6 . This is not correct since T n oc E. The reason 
is that as E ->■ 0, (dk 2 /dx) a ->■ while (d 2 k 2 /dx 2 ) a = -2a is constant. 
This means that (6) is not valid and we cannot approximate solutions by 
Airy's functions. In order to correct this, we need to consider specifically the 
potential given by (4). Now, as E — > 0, k 2 — > —ax 2 for < x < b. Use a 
change of variable y = a x ^x, (1) becomes: 

iP yy -y 2 iP = 0. (22) 

This equation can be solved by so called parabolic cylinder functions 10 . Ac- 
tually, we may use these functions to solve (1) with k 2 given by (4) exactly 
and write the tunneling probability in closed form. This is out of the scope 
of this paper. So let us assume that we do not know these functions. We 
will see that we do not need to solve (22) exactly. Instead, we may study it 
by WKB method! First, note that (dy/dy) jy 2 = 1/y 2 <C 1 as y — > oo, i.e. 
the WKB condition (3) is satisfied asymptotically. So (22) has asymptotic 
solutions given by (2) with k 2 = —y 2 . Let us define two solutions of (22) by 
their asymptotic behaviors as y — > oo : 
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Then we may write the solution for < x < b as 

1> = [Bj'fa)Aj(y) - Aj'fa)Bj(y)] /2 , (24) 

where 77 = a 1//4 6, in order to match with an outgoing solution e ty ^ x — > 1 
for x > b. Note that we have used the fact that the Wronskian W = AjBj' — 
Aj'Bj = 2 and that E — >■ 0. Using similar steps as (19) to (20) and neglecting 
Aj'(^) as compared with Bj'(?y), the tunneling probability can be found as: 

To = 16£ /{a 1 / 2 [Aj'(0)Bj'(??)] 2 } . (25) 

This gives the correct E dependence since a and b are independent of E. To 
find Bj'(?7), we may use (23) as a first approximation: 

Bj'(y) « ^ 2 /2. (26 ) 

So, 

To = «v^;/>j'(o)] 2Tw - (27) 

This gives the dependence on E, a and E b since Aj'(0) is only a constant. 
In order to compare T with T n numerically, we need to know Aj'(0), which 
can be found by numerical integration of (22). I found 11 Aj'(0) w —0.9777. 
Using this, some values of T are shown in Table II. We can see that they are 
quite close, although there are more than 10% difference. The main error 
can be shown to be due to the approximation in (26). To see this, we note 
that the factor before the exponential function in (23) should actually be an 
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asymptotic series 10 . The series of Bj can be found by requiring cancellations 
between terms when it is put into (22): 

B ^-9( 1 + ^{ 1 + ^[ 1 + ^ (1 + -'# (28) 

This series can be evaluated numerically up to a term with smallest magni- 
tude. This brought a factor of 1.133 to T in Table II, e.g. T changed from 
5.97 x 1(T 10 to 6.76 x 10 -10 for E = 10 -8 . This is very close to T n which is 
6.78 x 1(T 10 . 

The dependence of a and E b in (27) were also be verified numerically by 
choosing different a and E b . In general, as a decreases or E b increases, i.e. 
T w — > 0, the difference between T and T n decreases. The fact that we can 
find out the dependence of the tunneling probability on E, a and E b by using 
WKB idea without solving (22) shows how powerful WKB method may be 
if applied correctly. 

VI. For E « E h 

For E w E b , (14) and (17) no longer give values close to T n as we can see from 
Table I. However, we still can calculate the tunneling probability by WKB 
method. Now, 6 — a is small, we may approximate (4) by k 2 — \ 3 (x — a) for 
x in a region near a, including b, where A 3 = 2(aE) 1 ^ 2 . Then for this region, 
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the solution is given by: 



= Tre^ 6 { [Bi'(0 - jVEBi(0 



Ai (z) + 



-Ai'(0 + jVEAm 



Bi(z) 



(29) 



where z = A(x — a), £ = A (6 — a), in order to match with an outgoing solution 



iVEx 



for x > b. Using standard connection formula (7), we may connect it 



to WKB type solutions (2) for x — > — oo. Then, by grouping the coefficient 
of the incoming terms, we found the tunneling probability: 

4Av / £/vr 



2 ' 



(30) 



[Bi(0v<E - AAi'(0] + [M(£)Ve + ABi'(0; ' 
The Airy's functions can be evaluated by power series expansions (21). Some 
values of T b are shown in Table I for E < E b and in Table III for E > E b . 
We see that T b gives a quite good approximation to T n for \E — E b \ < 0.7. 
The range of validity for this approximation is surprisingly large at first 
sight. However, if we remember that the potential is proportional to x 2 so 
that although \E — E b \ is not small, \a — b\ may be small enough for the 
approximation to work. 

VII. For E > E b 

For the case that E larger than E b , T w — 1 while the true "tunneling prob- 
ability" (it may be better to call it transmission coefficient now) may differ 
from 1 quite a lot. However, WKB method still gives a good approximation. 
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We now assume that WKB condition (3) is satisfied for x < b all the way 
to — oo. Then, we only need to match the outgoing solution e % ^ x for x > b 
with WKB type solutions (2) for x < b. Find the coefficient of the incoming 
term. Then the transmission coefficient can be found 



AJE(E - E b ) 

= r / ; 1 , . (si) 

\/E — E b + \[E 

Some values of Too are shown in Table III. Since both and T n — >■ 1 as E — > 
oo. It is more appropriate to compare reflection coefficients R^ = 1 — 
and R n = 1 — T n . We see that i?oo gives a very good approximation to R n 
for large E. For E close to E b) fails and we need to use T b instead. 

VIII. Conclusions 

From the above discussions, we see that the standard WKB tunneling factor 
T w fails badly for E — > 0, E ~ E b and may have error up to factor of two in 
between. However, WKB approximation not necessarily fail provided that 
we use appropriate connection relations. For the example that k 2 given by 
(4), WKB approximation works for almost entire energy range from E — 
to E — > oo if we use different connection relations for different ranges of 
energy. Our conclusion is that when we use the standard WKB formula T w , 
we need to be very careful. If higher accuracy is desired, we need to consider 
connection relations case by case. 
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Appendix: Numerical calculation of T n 



Although a numerical method to find T n was described in ref. 4, a different 
method was used in this paper. Let 

^(&) = e ^ h 
Mb) = iVEe^ b , 

so that ip connects to an outgoing solution e l ^ x for x > b. Use Runga-Kutta 
method to integrate (1), with k 2 given by (4), from x = b back to x — 0. 
Match with incoming and reflected waves for x < and find the coefficient 
of the incoming term. Then the transmission coefficient can be found by 

2 



T„ = 4 



15 



Table I. 



Comparison of tunneling probabilities, T w ,Ti, T 2 and T b calculated by WKB 
method using different approximations with T n calculated numerically, for 
energy E in the middle range. Parameters a = 0.040965, E b = 1.2776 were 
chosen so that we may compare with Table I of ref. 4. 



E 


T 




T 2 


T b 


T 





0.00181 













0.1 


0.00611 


0.00656 


0.00654 


0.0342 


0.00987 


0.2 


0.0144 


0.0210 


0.0208 


0.0543 


0.0285 


0.3 


0.0297 


0.0504 


0.0491 


0.0896 


0.0609 


0.4 


0.0555 


0.103 


0.0979 


0.142 


0.112 


0.5 


0.0962 


0.188 


0.171 


0.214 


0.187 


0.6 


0.156 


0.312 


0.269 


0.301 


0.283 


0.7 


0.240 


0.479 


0.382 


0.397 


0.391 


0.8 


0.351 


0.678 


0.495 


0.494 


0.500 


0.9 


0.486 


0.886 


0.590 


0.584 


0.597 


1.0 


0.640 


1.06 


0.647 


0.663 


0.678 


1.1 


0.799 


1.11 


0.646 


0.729 


0.741 


1.2 


0.938 


0.896 


0.537 


0.783 


0.791 


E b 


1 








0.817 


0.821 
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Table II. 



Comparison of tunneling probabilities, T w , T 2 and T , calculated by WKB 
method using different approximations with T n calculated numerically, for 
energy E {a = 0.040965, E b = 1.2776). 



E 




T 




T 2 




T 




T 




10" 


-8 


1.81 x 10- 


-3 


6.42 x 10- 


-7 


5.97 x 10" 


-10 


6.78 x 10- 


-10 


10" 


-7 


1.81 x 10- 


-3 


2.03 x 10- 


-6 


5.97 x 10" 


-9 


6.78 x 10" 


-9 


io- 


-6 


1.81 x 10- 


-3 


6.42 x 10" 


-6 


5.97 x 10- 


-8 


6.78 x 10" 


-8 


10" 


-5 


1.81 x 10- 


-3 


2.03 x 10- 


-5 


5.97 x 10- 


-7 


6.78 x 10- 


-7 


io- 


-4 


1.82 x 10- 


-3 


6.44 x 10- 


-5 


5.99 x 10" 


-6 


6.78 x 10- 


-6 


10" 


-3 


1.86 x 10" 


-3 


2.08 x 10- 


-4 


6.11 x 10- 


-5 


6.80 x 10- 


-5 


10" 


-2 


2.17 x 10- 


-3 


7.64 x 10- 


-4 


7.11 x 10- 


-4 


7.04 x 10- 


-4 


10" 


-1 


6.11 x 10- 


-3 


6.54 x 10- 


-3 


1.93 x 10" 


-2 


9.87 x 10" 


-3 
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Table III. 



Comparison of tunneling probabilities, T& and T^, calculated by WKB 
method using different approximations with T n calculated numerically, for 
energy E > E b . = 1 — T^and R n = 1 — T n are reflection coefficients 
(a = 0.040965, £ fe = 1.2776). 



E 




T 


T 




Rn 


1.3 


0.826 


0.410 


0.829 


0.509 


0.171 


1.6 


0.907 


0.855 


0.901 


0.145 


0.0986 


2.0 


0.956 


0.938 


0.947 


0.0621 


0.0528 


2.6 


0.985 


0.972 


0.974 


0.0280 


0.0262 


3.0 


1.01 


0.981 


0.982 


0.0190 


0.0184 


4.0 




0.991 


0.991 


0.00920 


0.00901 


5.0 




0.995 


0.995 


0.00542 


0.00539 


6.0 




0.996 


0.996 


0.00357 


0.00356 


7.0 




0.997 


0.997 


0.00253 


0.00252 
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